NASA/TM— 2005-213821 


AIAA-2005-2941 



Effect of Free Jet on Refraction and Noise 


Abbas Khavaran 

QSS Group, Inc., Cleveland, Ohio 

Nicholas J. Georgiadis, James E. Bridges, and Vance F. Dippold III 
Glenn Research Center, Cleveland, Ohio 


July 2005 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/TM— 2005-213821 


AIAA-2005-2941 



Effect of Free Jet on Refraction and Noise 


Abbas Khavaran 

QSS Group, Inc., Cleveland, Ohio 

Nicholas J. Georgiadis, James E. Bridges, and Vance F. Dippold III 
Glenn Research Center, Cleveland, Ohio 


Prepared for the 

11th Aeroacoustics Conference 

cosponsored by the American Institute of Aeronautics and Astronautics 
and the Confederation of European Aerospace Societies 
Monterey, California, May 23-25, 2005 


National Aeronautics and 
Space Administration 


Glenn Research Center 


July 2005 



Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / /gltrs. grc.nasa.gov 


Effect of Free Jet on Refraction and Noise 


Abbas Khavaran 
QSS Group, Inc. 

Cleveland, Ohio 44135 

Nicholas J. Georgiadis, James E. Bridges, and Vance F. Dippold III 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


This article investigates the role of a free jet on the sound radiated from a jet. In particular, the role of an infinite 
wind tunnel, which simulates the forward flight condition, is compared to that of a finite wind tunnel. The second 
configuration is usually used in experiments, where the microphones are located in a static ambient medium far 
outside the free jet. To study the effect of the free jet on noise, both propagation and source strength need to be 
addressed. In this work, the exact Green’s function in a locally parallel flow is derived for a simulated flight case. 
Numerical examples are presented that show a reduction in the magnitude of the Green’s function in the aft arc and 
an increase in the forward arc for the simulated flight condition. The effect of finite wind tunnel on refraction is 
sensitive to the source location and is most pronounced in the aft arc. A Reynolds-averaged Navier-Stokes solution 
(RANS) yields the required mean flow and turbulence scales that are used in the jet mixing noise spectrum 
calculations. In addition to the sound/flow interaction, the separate effect of source strength and elongation of the 
noise-generating region of the jet in a forward flight is studied. Comparisons are made with experiments for the 
static and finite tunnel cases. Finally, the standard free-jet shear corrections that convert the finite wind tunnel 
measurements to an ideal wind tunnel arrangement are evaluated. 


I. Introduction 

Forward flight is generally believed to reduce the jet noise emission due to the reduced shear. Experimentally, 
simulation of the flight and its effect on jet noise is carried out in either an open wind tunnel (OWT) or an ideal wind 
tunnel (IWT). In the first configuration, the primary jet is surrounded by a co-flow (free-jet) with Mach number M* 
that simulates the flight Mach number, while the microphones are positioned outside the free stream in the stationary 
atmosphere. In the second arrangement, also referred to as an infinite wind tunnel, the microphones are positioned 
inside the co-flow to avoid the ambiguities associated with the refraction correction. In order for the microphones to 
be in the far field of the acoustic sources, the tunnel flow needs to extend far enough in the span-wise direction. 
Obviously, in either case, the measurements are carried out with fixed microphones. 

In practice, the OWT measurements are corrected for the secondary shear layer refraction, using semi-empirical 
corrections equations. The usual practice (ref. 1) is to replace the free jet shear layer with a vortex-sheet and to 
consider the source along the jet center line-and use a three-step calibration procedure for angle, distance, and sound 
amplitude. This procedure takes advantage of the high frequency geometric acoustics arguments, i.e., conservation 
of acoustic energy along a ray tube, and converts the noise measured outside the free jet to the levels that would be 
measured in an IWT. These simplifications, among other factors, ignore the difference between the on-axis and off- 
axis source location (ref. 2) as well as the additional noise generated by the free-jet shear layer. 

Having completed the IWT simulation, the actual flight conditions may be addressed correctly only when the 
observer is set in motion relative to the stationary nozzle, at a velocity that equals the IWT velocity. The frequency 
and the spectral density measured by such a moving microphone are related to those observed by the stationary 
microphone. This final step recreates the relative motion of an airplane with respect to a stationary observer on the 
ground. 

Most semi-empirical predictions methods (refs. 3 and 4) exploit the usual Lighthill’s V s scaling law and write 
the overall sound pressure level in flight as V" id V" , where V re i is the jet exit velocity relative to the ambient. The 
flight exponent factor m is allowed to change from 5 to 7.2 depending on the angle, and factor n is close to 3.0. 
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Michalke and Michel (refs. 5 and 6) used a convected form of Lighthill’s equation and proposed scaling laws 
that calculate the reduction level of the sideline sound intensity from static to flight. Tanna and Morris (ref. 7) 
examined the noise reduction in the forward flight from unheated jets, and suggested that a relative flight exponent 
factor of 5.5 shows reasonable agreement with the 90° measurements. Tam et al. (ref. 8) studied the influence of 
forward flight on the sources of fine-scale turbulence mixing noise in an OWT simulation. They calculated the mean 
flow using specified velocity and temperature profiles at the jet exit that simulated an OWT case, and demonstrated 
a reduction in the jet turbulent kinetic energy and its radiated noise along the sideline. 

The reduced shear due to the forward flight velocity affects not only the source intensity and its spectral 
distribution (ref. 9), but also the sound propagation. The purpose of the present work is to present a prediction 
method that accounts for source modification as well as refraction of sound due to flight and to provide a measure of 
comparison between the IWT and OWT simulations relative to the static case. The Green’s function (GF) for a static 
ambient medium (refs. 10 and 11) is now extended to the IWT condition via the adjoint method-and is calculated 
exactly in an axisymmetric jet. The subsequent spectral predictions show that noise reduction in flight is not uniform 
throughout the spectrum at all angles. There is less high frequency reduction at small aft angles. Examination of the 
source spectral density demonstrates that flight tends to stretch the noise-producing region of the jet plume and 
reduce its peak intensity level. These results and in particular the mean flow effect on propagation from static to 
OWT and IWT cases will be explored and compared with the usual shear layer refraction correction. 

The paper is organized in the following order. The Green’s function formulation in an IWT simulation is 
described in section 2 and is utilized in jet noise calculations described in section 3. Sample numerical examples 
using a closed-form mean velocity profile are shown in section 4. Details of the GF in numerically computed jets are 
addressed is section 5, and the GF is mapped for the entire jet at selective observer angles and frequencies in static 
as well as in IWT and OWT simulations. The influence of forward flight on both the source spectral density and the 
stretching of the noise producing regions of a jet is studied in section 6. Section 7 examines the far- field noise in 
various simulations and compares those results with the OWT data as well as with the refraction-corrected data 
using a guideline commonly practiced in industry. A transformation from the IWT simulation to the actual flight 
condition is described in section 8. Section 9 provides the final summary and conclusions. 

II. Formulation 

The adjoint Lilley equation governing a harmonic point source with frequency ftiis 

L a [^ a (x,x°,co)e- ica ] = cy M S(x-x°) (1) 

where L a denotes the adjoint operator in a locally parallel flow (ref. 11), and S fl is the adjoint Green’s function 
which relates to the direct Green’s functional by switching the source and observer locations 


d(x° ,x s ,a>) = d a (x s ,x° ,a>). 

Consider the case when the main jet is surrounded by an infinite free-stream. Outside the jet, we should have 

U = U X , p = p m , P = P^ as r — > co . 

In the ambient medium, equation (1) becomes 


a 


a 


-a CO + U K , y+d(ico + U—)V l \% a = ci,d(x -x°). 

I UJi i yJJv J | 

The solution to equation (4) that satisfies the outgoing wave conditions is 

S = exp{/'— \x-x °\) , a> = — , 

A7tia>\x-x°\ c 1+M cos<9 


(2) 


( 3 ) 


( 4 ) 


( 5 ) 
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Where M x = UJc x is the ambient Mach number and the polar angle 0 is measured from down stream jet axis, 
COS# = (x° — X t ) / \x° — x|. Wave number k is now defined as 

k=mic x (6) 

Upon introducing the spherical coordinates ( R , 0. (p), with R measured from the jet exit centerline and the azimuthal 
angle ^measured in a span-wise plane, we write equation (5) as 

= l —^z — exp{-z'£(xj cos# + x 2 sin#- R)}, r (7) 

4 nc K k R 

The above plane wave expression may be written as a superposition of cylindrical functions 

g-ik{x x cos6-R) oo 

e fl =— ; X (-0" e m J m (kr sin #) cosm((Z> -cp 0 ), r -> go (8) 

4 ncJiR m= ci 


In the free stream r > R 0 , U=U^ and the scattered part of the Green’s function is written in terms of wave number k 
as 

-i k (xyCOS0-R) oo 

^ =~ A f —Y J h m (r,k,0)cosm((p-(p o ), r>R o (9) 

4 ncJtR m=0 

The adjoint GF in the far- field is a sum of the incident and scattered waves 

6 a= 6 a + 6 « ) 5 r ^ R o ( 10 ) 

r-» oo 

Substituting equation (10) into equation (4) and upon noticing that equation (8) already satisfies equation (4), then 
equation (9) should satisfy the homogeneous form of equation (4) 

=o (ii) 

Equation (11) has a solution as the mth-order Hankel function \k r sin #) 

-ik(x x c,o$G-R) oo 

= — ~^Z A m H m(krsrnO ) cos m{(p - (p° ) , r > R 0 (12) 

4 nc^kR m=0 

Using equation (8) and equation (12) into equations (10), the outer solution is now matched with the inner solution 

g-i k (xj o,osO-R) oo 

& a (x s ,x° ,co) = &(x° ,x s ,co) = — Z f m (r,k,0)cosrn((p-<p o ), r < R o (13) 

4 ncJkR m=0 

at the jet boundary r = R a , i.e., the two solutions and their radial derivatives are set equal. 

After substituting the inner solution equation (13) into the homogeneous form of equation (1), the function f m is 
found to satisfy the compressible Rayleigh equation 
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fm + 


— - — - 4 cos$- 
r p 


M' 


( 1 - McosO) 


fm + 


0) ~ „ 2 m CO 2 n 3 cos# .p' ~ 1 ~ ~ . 

— r (l- McosQy t — cos' G + ~ (— M' — M -M ) 

c r c; (1-Mcos#) p r 


( 14 ) 


\L = 0 , 


where 


M = M - M m , and M = 


U(r ) 


(15) 


In equation (14), prime denotes a radial derivative. Notice that for a given observer frequency to, the above adjoint 
equation is solved in a frame of reference moving with U x and CO as defined in equation (5) is the frequency 
observed in this reference frame. As r— »oo, M as well as M' and p' approach zero and equation (14) reduces to 
Helmholtz equation (11), thus stratifying the matching conditions at the jet boundary. 

The numerical solution to (14) starts with /„, = a m r m as r— >0, and continues to the jet boundary r = R 0 where the 
matching conditions are applied to determine A m and a m 

a , , y . £ iih iAiQCM 

m CsmfffJRJN"«(z)-d/JRJ/c/rH< m ‘>( Z ) (16) 

/ = sin ft - 1 2 - o I , m = 0,1,2,... 

In equation (16), prime as appearing on the Bessel and Hankel functions points to a derivative with respect to the 
argument. 


III. Application to Jet Noise 

The above Green’s function is now extended to a moving type singularity of interest in jet noise. The governing 
equation is Lilley’s third-order wave equation. This equation is linearized about a unidirectional transversely sheared 
mean flow (ref. 12) 

p„= constant, U, = SnU(x 2 *3) T 0 =T 0 (x 2 x 3) 


Dt dx ! ax, ax, 

Ct =^(i+^'Kw y +(r-i) h '-^- 

ox, ox i 

J 1 


(17) 


(18a) 


L is linearized Lilley operator, Q denotes the equivalent source component in a Lilley-type acoustic analogy, and the 
convective derivative is D I Dt = d / dt + Ud / dx x . X = (x x ,X 2 ,X 3 ) are the Cartesian coordinates and x t is in the 

stream-wise direction, t is time, tt'=(p/ p o Y r ~ 1, p is pressure and p 0 is the mean pressure, (i.e., 
p = p o + p' ), T„ is the mean temperature, m, is the fluctuating velocity component u-, = v— 5,i U. h' denotes enthalpy 

fluctuation h'= (/ — l)(c~)', which is related to fluctuations in sound speed (c 2 )'= y^KT— T g ) . Since pressure 
fluctuations p' are small relative to mean pressure p u . the dependent variable is approximated as 7l' = p'/( V p o ) . 
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In isothermal jets, the source term associated with enthalpy fluctuations is neglected, and equation (18a) is 
approximated as 


^ dU'Uj 
8X: 


(18b) 


The two source terms on the right hand side of equation (17) are usually referred to as self- and shear noise terms. 


Consider a convecting harmonic point source with a source frequency 0 ) s and convection velocity i U c . The GF 
associated with the first source term in equation (17) is 


L{Ge 


Dt 


2 e imt 8{x l -U c t)S{x t -x s t )} . 


(19) 


Subscript t denotes a transverse location. The above GF is derived from the stationary simple-source GF (eq. 13), 
and the convolution integral. Following the steps described in (ref. 10) it is shown that 


G(x,x s ,a>) 


-i (1- M s cos 6) 
4 xR (1 - M c cos #) 


e ' kR Y,fnX rS > k , #) cos m(<p - (p s ). 

m = 0 


(20) 


Wave number k was defined in equation (6) and 


M s = U (r s ) / c x , M c = U C / c x , M C = M C -M X . 

In deriving equation (20), the requirement that the time harmonic factor exp {-icot) should cancel out from both sides 
of the equation results in the following definition for the Doppler factor 


a> s ( 1+ M x cos#) 

CO = 

1 - M cos# 


(21) 


Since M c —> 0 as r— >oo, then ©— »co s in the infinite free stream. 


In axisymmetric jets, the source intensity is constant along a ring in a span-wise plane at (A',' ,r' ). A ring source 
directivity factor is defined below that is used in conjunction with the source/GF volume integration 


D 2 (x,x s ,co) = 


(4 tzR) 2 
2k 


j\G(x,x s ,a))\ dtp s . 

-71 


(22) 


IV. Numerical Example 

To demonstrate the effect of the free-stream flow on the Green’s function, consider a jet profile described as 


Pir) _ x 

Poo 

^± = M m +(Mj- M x )sech 2 (2r/Dj). 


(23a) 
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The above mean velocity profile is shown in figure 1 using a jet Mach number of Mj = 0.90. Two profiles using the 
static Mrc = 0.0 and the ideal wind tunnel conditions M„ = 0.18 are examined. 

A third profile is also shown that is obtained by subtracting 

^l = M 00 {l-sech 2 (2r/ J D J -4)}, r>2D, (23b) 

Coo 

from equation (23a). This profile, seen as dash-dot, decays from a free-stream value of M x = 0.18 to 0.0 at 
r > ID j, thus simulating an open wind tunnel condition. 

The directivity factor D as defined in equation (22) was calculated numerically for the three mean flow profiles 
shown in figure 1. Figures 2 to 5 show the effect of the mean flow on a stationary ring source ( U c = 0.0), at a 
Strouhal number St = coDJ{2nUj) of 0.50. The static case M, = 0.0, as well as the IWT and OWT simulation is 
shown at a flight Mach number of M x = 0.18. Figure 2 shows that in the vicinity of the peak directivity angle, and 
prior to the zone of relative silence, there is stronger refraction and beaming of sound arising from a two-steps layer 
(i.e., OWT) compared to a single shear layer (i.e., IWT). As the source moves away from the centerline, the OWT 
shear layers could surpass even the static case in its sound beaming in the direction of the peak directivity angle 
(figs. 3 and 4). In the forward arc, for the most part, there is little difference between these two profiles - and both 
simulations appear to amplify the sound relative to the static case at very large forward angles. 

In general, the effect of the infinite wind tunnel, compared to the static conditions, may be surmised from 
figures 2 and 3, as a reduction of noise in the aft arc and an increase in the utmost forward arc. Thus, the infinite free 
stream plays opposing roles in the forward and aft arcs. 

It is noted that when calculating the actual jet noise spectrum, source strength needs to be convolved with an 
appropriate GF. Turbulence intensity, and hence the source strength, are affected by a reduction in the mean shear 
due to the presence of the free jet. This should naturally lower the jet noise level at 90° even though the isothermal 
GF appears unaffected as discussed above. 

V. Mapping of the Green’s Function 

To further illustrate the mean flow effect on the radiated sound in a jet, it is of interest to study the GF in more 
detail as a function of the source location and jet profile. For this purpose, we consider a numerically computed jet at 
Mj = 0.98, T r = 1.0 (i.e., Uj / c w = 0.90), and select secondary flow streams that would simulate both the IWT and 
OWT conditions. Table 1 shows the secondary Mach, or for simulations that follow. 


TABLE 1.— MACH 0.98 UNHEATED JET 


IWT Mach no. 

OWT Mach no. 

0.00 


0.10 

0.10 

0.20 

0.20 

0.28 

0.28 


Mean flow predictions (RANS) were generated using the WIND flow solver (refs. 13 and 14) available at the 
NASA Glenn Research Center. The inflow conditions are specified at five diameters upstream of the jet exit plane, 
and in terms of plenum total pressure and temperature. In order to establish a given Mach number for the outer flow, 
the static pressure and the desired free-stream Mach number are used to calculate the inflow total pressure. 
Similarly, the inflow total temperature is defined from the static temperature and the free-stream Mach. While in an 
IWT simulation, these two conditions are set as constant throughout inflow boundary, an OWT simulation requires a 
gradual reduction in the inflow total pressure and temperature to their respective static ambient values in an interval, 
say 5.5 D, to IDj. This method is advantageous over an alternative that would specify inflow velocity and 
temperature profiles-as it closely simulates the actual test conditions, allows for the secondary entrainment, and lets 
the required mean velocity and temperature profiles in an OWT simulation to develop ahead of the nozzle exit. The 
presence of the ambient flow, in either simulation, elongates the length of the potential core and reduces the peak 
turbulent kinetic energy (see figs. 6 and 7). 
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The RANS grid extends 40 diameters downstream of the nozzle exit plane, and 20 diameters in the span-wise 
direction. The mean axial velocity profiles (fig. 8) are shown for the static as well as the M* = 0.20 (IWT and OWT) 
simulations. Similar illustrations are presented for the turbulent kinetic energy profiles (fig. 9). At a radial distance 
of 10 Dj (not shown in fig. 9), the turbulent kinetic energy level decays throughout the jet to 0.0, 0.80, and 
2.5 percent for the static, IWT, and OWT simulations respectively. Between the IWT and OWT simulations the peak 
turbulence levels are very close, consequently the additional free-stream turbulence (larger scales) tend to produce 
relatively more low frequency noise in the OWT case as will be shown in section 7. 

The ring directivity factor D in equation (22) was calculated at selective inlet angles of (60°, 90°, 120°, and 
150°) throughout the computational domain starting from the jet exit plane. Using a locally parallel flow model, the 
mean temperature T(x,r) and axial velocity U(x,r) were used in the equation (14). The number of modes required to 
achieve a converged solution in equation (20) would increase with frequency. Results are presented for a stationary 
source (M c = 0.0) at three Strouhal numbers of St = (0.10, 0.50, 3.0), and using maximum mode numbers of 
to = (7, 10, 20), respectively. Three simulations are shown: (1) static M m = 0; (2) IWT with M w = 0.20; and (3) OWT 
with M x = 0.20. 

The Green’s function D as defined in equation (22) is mapped as a function of source location at St = 0.10 
(fig. 10) at four angles in each simulation. Identical color legends are shown at equal angles in order to highlight the 
differences in the three simulations (seen side-by-side). The static GF could have an amplifying effect on the sound 
radiated to 150° if volume sources are located near the centerline at 4 to 7 diameters from the jet exit. Fortunately 
there is very little turbulent kinetic energy present in this vicinity, which coincides with the jet potential core 
(see fig. 7). Between the IWT and OWT simulations, the latter demonstrates stronger low frequency directionality at 
120° and 150° observer angles. 

Figures 11 and 12 exhibit similar comparisons at St = 0.50 and 3.0, respectively. At shallow aft angles, the 
OWT flow achieves a more intense high frequency refraction compared to its IWT counterpart (see fig. 12). 
Subsequent spectrum predictions confirm this observation. 

VI. Source Spectral Density 

Flight reduces the peak turbulent kinetic energy and stretches the noise-producing region of a jet. It also affects 
the time-and length-scales of a noise generating turbulent eddy. The far-field spectra from source terms of 
equation (17) have been studied under the static condition (ref. 10). Those results are readily extended to a simulated 
flight case (app. A) using the GF definition provided in equation (20). The far-field sound at each observer angle and 
frequency of interest is now calculated from a volume integration of equations (Al) and (A2) over the entire jet 
volume. 

Sound spectral density, when calculated per unit axial slice at each stream-wise grid coordinate, provides 
valuable insight into the source strength and its spectral distribution as affected by flight. For simplicity, consider a 
90° far-field observer at RJD t = 100. The sound spectral density (eq. Al) was calculated in a simulated flight at an 
observer Strouhal number range of 0.016 to 18.0. 

Figure 13 shows the source spectral distribution in the first 25 diameters of the jet, using free jet Mach numbers 
of 0.0, 0.20, and 0.28 in an IWT simulation. Contour legends extend 20dB from the peak in each simulation. High 
intensity noise producing sources move downstream in excess of 3 diameters as M x in increased from 0.0 to 0.28. 

An OWT simulation works similarly to reduce source intensity and shifts its maximum location. For example, 
figure 14 shows that there is very little difference in source distribution between the IWT and OWT simulations as 
viewed by an observer at 90°. Any spectral difference between the two simulations is primarily linked to the Green’s 
function. 


VII. Far-field Sound 

Spectrum predictions are presented for Mach 0.98 unheated jets issuing from a 5.08 cm (2.0 in.) convergent 
nozzle on an arc at 100 Dj. The IWT simulations using ambient conditions of table 1 are shown in figure 15. In the 
mid-angle range, there are 2.5- to 3 dB uniform reductions in sound level across the spectrum from one simulation to 
the next. At shallow angles, however, these reductions appear to diminish at high frequency (HF) as the tunnel Mach 
number is increased. Measurements shown in figure 15 are acquired at the Small Hot Jet Acoustics Rig (SHJAR) at 
the NASA Glenn Research Center, and are relevant to the static condition only (ref. 15). 

As pointed out earlier, flight simulation measurements are usually carried out in an OWT flow, and are 
subsequently processed for the secondary shear layer refraction using semi-empirical correction relations. Figure 16 
shows a comparison between the IWT and OWT far-field noise predictions at = 0.10. Two sets of data are also 
shown; (a) the actual measurements in an OWT simulation (ref. 16); (b) the IWT data that are obtained from the 
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former set after utilizing the usual refraction correction recipe (ref. 1). All data are converted to a narrow band and 
are presented as lossless. 

The differences in spectral levels between the two simulations become increasingly visible at aft angles. Both 
data and prediction indicate that at low frequency the OWT flow is noisier compared to its corresponding IWT flow. 
There is less agreement in the HF range however. Predictions suggest less HF noise at shallow angles in an OWT 
flow (relative to an IWT flow) due to a stronger F1F refraction. The correction procedure remains rather insensitive 
to this observation. Similar conclusions can be drawn at M x = 0.20 as seen in figure 17. 

VIII. Influence of Flight 

A transformation to the actual flight conditions can be carried out if the microphone is set in motion at M x 
within the IWT simulation. Such a microphone measures a frequency CO ' that relates to co as measured by a 
stationary microphone according to CO' = 0)(\+M x COS#) 1 . Since the over all power level measured in either 
frame remains the same, the sound spectral intensity should change as G)' /' (co' ) = CO I (co) . This results in a 
dynamic factor that relates the spectral densities according to 

p' 2 (co') = (1 + M x COS 6) p 2 {co) . (24) 

Flere prime refers to a moving microphone within an IWT - that would also represent a stationary observer 
listening to a fly-by airplane. Note that the spectral density on the right-hand-side of equation (24) is not that of a 
static condition, rather it is due to a stationary jet engine in an infinite wind tunnel simulation. Measurement of the 
actual jet mixing noise in flight is further complicated due to the presence of other types of sources, such as 
boundary layer noise produced by the flow around the engine nacelle and the external surfaces, and possible 
interaction between different types of sources. This might explain the commonly observed differences between the 
flight simulation experiments and flyover tests. Most semi-empirical relations proposed in the literature relate the 
static engine noise (i.e., = 0) to the actual flight noise. These models suggest a wide range of flight dynamic 

factors (refs. 3, 5, 7, and 17) that account for the jet stretching effect (fig. 13) and the reduction in source strength 
due to flight in one form or another. 

Figure 18 shows the influence of the dynamic factor as defined in equation (24) on the sound spectral level and 
its frequency using a flyover Mach number of = 0.20. A comparison is also made with results shown earlier for a 
fixed microphone in an IWT simulation. 

IX. Summary 

In this work, the effect of forward flight on the high-speed jet noise was investigated within the framework of a 
Lilley-type acoustic analogy. The influence of flight on the source spectral intensity and sound propagation was 
addressed in a systematic order. Both ideal wind tunnel and open wind tunnel simulations were studied using 
carefully specified RANS solutions that emulated the actual test conditions. It was argued that an IWT simulation 
more closely resembles the flight conditions, and reduces the spectral level in flight almost uniformly at sideline 
angles. At shallow aft angle, however, the high frequency noise reduction diminishes gradually as the flight Mach 
number is increased. 

Theoretical predictions demonstrate that the difference between the two simulations becomes increasingly 
visible at smaller aft angles. This is evidenced by an increase in the low frequency noise and a decrease in high 
frequency noise in an OWT simulation as compared to its IWT counterpart. The commonly practiced refraction 
correction procedure appears to remain insensitive to this predicted F1F noise difference between the two 
simulations. 

Flight, moreover, moves and stretches the noise producing regions of a jet in the direction of the jet flow. This 
could have an impact on the performance of noise silencer kits in the static versus flight conditions. Transformation 
from an IWT simulation (with static microphones) to the actual flight resulted is a flight dynamic factor that 
amplified the noise in the aft angles and additionally modified the observer frequency due to the flight Doppler 
factor. 


Appendix A. Far-field Spectral Density 

Expressions for the far-field noise due to the source terms of equation (17) were provided in reference 10 
assuming an isotropic turbulence, and using exponential functions to model both spatial and temporal dependencies 
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of a two-point space-time correlation. The GF for a simulated flight as given in equation (20) is used to extend the 
previous static results to a simulated flight case 
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The 4 th order axial correlation coefficient at source point V is 

I nu (y,a> s ) = A,y n 


\ + {co s z 0 / 2)' 


-N(k£) 


(A3) 


of is the source frequency which is related to the observer frequency according to equation (21), and 
l = CqK 3/2 Is, X 0 = C Z K / £ , where k and s are the turbulent kinetic energy is its dissipation rate. Parameters 
(c f ,C T ,A m ) denote three calibration constants, and the source convection velocity is defined as 
U c = 0.50U(x,r)+0.25Uj. 


Factor N(k€) in (A3) represents a non-compactness factor that is determined from the spatial function of the 
correlation. Its argument kf = ( C t / C r )(<yz" 0 )(/f 0 5 / C^) , as discussed in (ref. 10), becomes self-similar when 

divided by the local acoustic Mach number U(r s )/c x . This modified argument is used in noise predictions. In the 
presence of a free jet with velocity Uf ree , (i.e., IWT or OWT simulations), U(r') should be measured with respect to 
the free-stream velocity as U{r s )-Uf ree . 
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Figure 1. — Mean velocity profile with = 0.0 (solid line); 

Mo, = 0.18 ideal wind tunnel (dashed line); open wind tunnel 
(dash-dot). 



Figure 2. — Effect of the free jet on refraction. Green’s Figure 3. — As figure 2 but at r s /Dj= 0.50. 

function D{x,x s ,co) at observer Strouhal number 
of St = 0.50 in an isothermal Mj = 0.90 jet. 

Source at r % tDj= 0.0 with: M x = 0.0, (solid line); 

M„o = 0.18 IWT, (dashed line); OWT (dash-dot). 
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Figure 4. — As figure 2 but at r s /Dj= 1.0. 


Figure 5. — As figure 2 but at rVDj = 2.0. 
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Figure 6. — Centerline velocity decay in a simulated flight, Mach 
0.98 unheated jet. Solid lines (IWT), from the top: = 0.28, 

0.20, 0.10, 0.0. Dashed lines (OWT), from the top: M m = 0.28, 
0.20, 0.10 



Figure 7. — Centerline turbulent kinetic energy in a simulated 
flight, Mach 0.98 unheated jet. Solid lines (IWT), from the 
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u/u, u/u, u/u 


bottom: M* = 0.28, 0.20, 0.10, 0.0. Dashed lines (OWT), 
from the bottom: M m = 0.28, 0.20, 0.10. 
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Figure 8. — Mean velocity profiles in a simulated 
flight, Mach 0.98 unheated jet. (a) M,,. = 0.0; 
(b) M» = 0.20 (IWT); (c) M x = 0.20 (OWT); at 
axial locations - starting from the top: 
x/Dj = 2, 5, 7, 10, 15,20. 


Figure 9. — Turbulent kinetic energy profiles in a 
simulated flight, Mach 0.98 unheated jet. 

(a) Moo = 0.0; (b) M» = 0.20 (IWT); 

(c) Moo = 0.20 (OWT). 
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Figure 10. — Green’s function mapped for a stationary monopole-type ring source in a Mach 0.98 unheated jet at 
St = 0.10. Inlet angles: from top 60°, 90°, 120° and 150°. Free jet simulations: from left M 00 = 0.0; M K = 0.20 
IWT; and M., = 0.20 OWT. 



Figure 11. — As figure 10 but at St = 0.50. 



Figure 12. — As figure 10 but at St = 3.0. 
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gure 13. — Source spectral density per unit axial slice in a Mach 0.98 unheated jet at 90°, 
and at R/Dj= 100.0. (a) M x = 0.0; (b) IWT, M m = 0.20; (c) IWT, M w = 0.28. 
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Figure 15. — Spectra in a Mach 0.98 unheated jet at indicated inlet angles and at R/Dj= 100. Ideal 
Wind Tunnel simulation atM x = 0.0, 0.10, 0.20, 0.28. Dashed line (prediction); solid line 
(static data-ref. 15). 
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Figure 16. — Spectra in a Mach 0.98 unheated jet at M x = 0.10 and at indicated inlet angles. Predictions: 
OWT (solid line); IWT (dashed line). Data: OWT (ref. 16), (open symbol); IWT correction 
(filled symbol). 
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Figure 17. — As figure 16 but at M m = 0.20. 



f Dj/Uj 


Figure 18. — Effect of airplane flyover on jet noise 
spectra in a Mach 0.98 unheated jet at M x = 0.2 
and at indicated inlet angles. Predictions: IWT 
(dashed line); IWT with moving microphone 
(solid line). 
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